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The NA62 experiment at CERN collected a large sample of ^ e+ v decays during a dedicated 
run in 2007, which allows a precise test of lepton universality by measurement of the helicity 
suppressed ratio Rk = r{K^ e+v)/r(/r+ M^v). The preliminary result of the analysis of 



a partial data sample of 51089 A'+ ^ e+v candidates is Rk 
with the Standard Model expectation. 



(2.500 ±0.016) X 10"^ consistent 
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Introduction 

The ratio of kaon leptonic decay rates Rk = T{K^ /T{K'^2) has been calculated with an 
excellent accuracy within the the Standard Model (SM) Rf^ = [mjm^f (^F^)^(l + 

5/?qed) = (2.477 ±0.001) x 10"^ where 5/?qed = (-3.78 ±0.04)% is a correction due to the in- 
ner bremsstrahlung (IB) K^jy process which is, unlike the structure dependent (SD) K^y process, by 
definition included into Rk- The factor {mf,/m^)^ accounts for the helicity suppression of the Ke2 
decay due to the V — A structure of the charged weak current, and enhances sensitivity to non-SM 
effects. In particular, enhancement of Rk by a few percent (relative) is possible in the MSSM with 
non-vanishing e — T lepton mixing [^, compatible to the presently known experimental constraints. 

The current world average composed of three 1970s measurements ^ and a recent KLOE 
result [§] is R^^ = (2.467 ±0.024) x 10"^. It has 1% relative precision and is compatible to the 
SM. The NA62 experiment at CERN collected data in 2007-08 aiming at an Rk measurement with 

0. 4% precision. The preliminary result obtained with a partial data sample is presented here. 

1. Beams, detector and data taking 

The NA48/2 beam line and setup were used; running conditions were optimized for the Ke2 
measurement in 2007 using the experience of earlier studies based on NA48/2 data sets [Q]. 

The beam line is capable of delivering simultaneous narrow momentum band and K~ 
beams; a central momentum 74 GeV/c was used in 2007. Momentum of the incoming kaon is not 
measured directly in every event; the beam average monitored with —>■ 3?!* decays is used to 
reconstruct K12 kinematics by missing mass Mmiss- A narrow momentum spectrum (Ap^^ / pK ~ 
2%) is used to minimize the corresponding contribution to the M^iss resolution. 

The K12 decay signature consists of a single reconstructed track. Since the incoming is 
not tracked, backgrounds induced by the beam halo have to be considered. The performance of 
the muon sweeping system results in lower background in sample (~ 1%) than in K^^ sample 
(~ 20%), therefore ~ 90% of the data were taken with the beam only, and small fractions were 
recorded with simultaneous beams and K beam only. The halo background is directly measurable 
using the samples of reconstructed K12 candidates of the sign not present in the beam. 

Among the subdetectors located downstream a vacuum decay volume, a magnetic spectrome- 
ter, a plastic scintillator hodoscope (HOD) and a liquid krypton electromagnetic calorimeter (LKr) 
are principal for the measurement. The spectrometer, used to detect charged products of kaon de- 
cays, is composed of four drift chambers (DCHs) and a dipole magnet. The HOD, used to produce 
fast trigger signals, consists of two planes of strip-shaped counters. The LKr, used for 7 detection 
and particle identification, is an almost homogeneous ionization chamber, 27Xo deep, segmented 
transversally into 13,248 cells (2x2 cm^ each), and with no longitudinal segmentation. A beam 
pipe traversing the centres of the detectors allows undecayed beam particles and muons from de- 
cays of beam pions to continue their path in vacuum. 

A minimum bias trigger configuration is employed, resulting in high efficiency with relatively 
low purity. The Kg2 trigger condition consists of coincidence of hits in the HOD planes (the so 
called Q\ signal) with 10 GeV LKr energy deposition. The trigger condition consists of the Qi 
signal alone downscaled by a factor of 150. Loose upper limits on DCH activity are also applied. 
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Most data taking took place during four months in 2007. Two weeks of data taking allocated 
in 2008 were used to collect special data samples for studies of systematic effects. The present 
analysis is based on ~ 40% of the 2007 data sample collected with the beam only. 

2. Analysis strategy and event selection 

Monte Carlo (MC) simulations are used to a limited extent only: 1) to evaluate a correction 
for the difference of Ke2 and K^2 geometric acceptances; 2) to simulate energetic bremsstrahlung 
by a muon, which is not directly accessible experimentally as discussed below. 

In order to compute geometrical acceptances, a detailed Geant3-based MC simulation is em- 
ployed. It includes full detector geometry and material description, stray magnetic fields, local 
inefficiencies, misalignment, detailed simulation of the beam line, and time variations of the above 
throughout the running period. The ^^2(7) processes are simulated in one-photon approximation |jl|]. 
Unlike the KLOE analysis the resummation of leading logarithms [|6|] is not included into the 
simulation for the preliminary result. 

The analysis strategy is based on counting the numbers of reconstructed Kei and K^2 candi- 
dates collected simultaneously, consequently the result does not rely on kaon flux measurement, 
and the systematic effects due to the detector efficiency cancel to first order. To take into account the 
momentum dependence of signal acceptance and background level, the measurement is performed 
independently in bins of reconstructed lepton momentum. The ratio Rk in each bin is computed as 

^ N{Ke2)-N^{K,2) UxAiK^2)xeiK^,2) 1 

D ' N{K^2) -A^b(^m2) ' feXAiK,2) X eiK,2) ' /lkt' 
where N{Ki2) are the numbers of selected Kg2 candidates {£ = e,pL), A^b(^£2) are numbers of back- 
ground events, fi are efficiencies of eIpL identification criteria, A{K(^2) are geometrical acceptances 
computed with MC, £{Ki2) are trigger efficiencies, /lki is the global efficiency of the LKr readout, 
and D = 150 is the downscaling factor of the K^2 trigger. 

Due to topological similarity of Ke2 and K^2 decays, a large part of the selection conditions is 
common for both decays: (1) exactly one reconstructed charged particle track; (2) its momentum 
15 GeV/c < p <65 GeV/c (the lower limit is due to the 10 GeV LKr energy deposit trigger 
requirement in Kg2 trigger); (3) extrapolated track impact points in DCH, LKr and HOD are within 
their geometrical acceptances; (4) no LKr energy deposition clusters with energy E > 2 GeV and 
not associated to the track to suppress background from other kaon decays; (5) distance between the 
charged track and the nominal kaon beam axis CDA < 1.5 cm, decay vertex longitudinal position 
within the nominal decay volume (the latter condition is optimized in each momentum bin). 

The following two principal selection criteria are different for the Kf,2 and K^2 decays. K£2 
kinematic identification is based on the reconstructed squared missing mass assuming the track to 
be an electron or a muon: M^^^^{£) = {Px—PtY, where Pa',^'^ {£ = e,/x) are the four-momenta of the 
kaon (average beam momentum assumed) and the lepton (electron or muon mass assumed). A cut 
|M^-j,j,(e)| < Mq is applied to select Ki,2 candidates, and \M'^^^^{}i)\ < Mq for ones, where Mq 
varies from 0.009 to 0.013 (GeV/c^)^ among track momentum bins depending on M^iss resolution. 
Particle identification is based on the ratio E/p of track energy deposit in the LKr to its momentum 
measured by the spectrometer. Particles with 0.95 < E/p < 1.1 {E/p < 0.85) are identified as 
electrons (muons). 
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Track momentum (GeV/c) Track momentum, GeV/c 



Figure 1: (a) Missing mass squared in electron hypothesis M^jj,j,(e) vs track momentum for reconstructed 
Ke2 and Ky^2 decays: kinematic separation of and K^2 decays is possible at low track momentum only, 
(b) Measured and simulated probability of muon identification as electron f (/i — > e) vs track momentum: 
data with the Pb wall, MC simulations with and without the Pb wall (signal region marked with arrows). 

3. Backgrounds 

K^2 decay is the main background source in the sample. Kinematic separation of Kg2 and 
K^2 decays is not achievable at high lepton momentum {p > 40 GeV/c), as shown in Fig. |l|a. The 
probability of muon identification as electron {E / p > 0.95 due to 'catastrophic' bremsstrahlung) 
is P{pi ^ ~ 3 X 10^^ in the NA62 experimental conditions, non-negligible compared to /?™ = 
2.477 X 10^^. Direct measurement of f (jU — > e) to ~ 10^^ precision is necessary for validation of a 
theoretical calculation of the bremsstrahlung cross-section [^] in the high y energy range, which is 
used to evaluate the K^2 background. Typical /i samples are affected by relatively large (~ 10^^) 
electron contamination due to /i ^ e decays in flight. To collect pure /i samples, a ~ IOXq thick 
lead (Pb) wall covering ~ 20% of the geometric acceptance was installed between the HOD planes 
during a part of the data taking. In the samples of tracks traversing the Pb and having E / p> 0.95, 
the electron component is suppressed to a level much below P{}i. —>■ e) due to energy loss in Pb. 

The momentum dependence of /'(/i — s- e) for muons traversing Pb was measured with a data 
sample collected during a special 20h muon run, and compared to the results of a dedicated Geant4- 
based MC simulation involving standard muon energy loss processes and bremsstrahlung according 
to The data/MC comparison (Fig. [I]b) shows excellent agreement in a wide momentum range 
within statistical errors, which validates the cross-section calculation at the corresponding precision 
level. The simulation shows that the Pb wall modifies e) via two principal mechanisms: 1) 

muon energy loss in the Pb by ionization decreasing P{}x e) and dominating at low momentum; 
2) bremsstrahlung in Pb increasing e) and dominating at high momentum. 

To estimate the K^2 background contamination, the kinematic suppression factor is computed 
with the standard setup simulation, while the validated simulation of muon interaction in the LKr 
is employed to account for P{pL — > e). The uncertainty of the background estimate comes from the 
limited size of the data sample used to validate the simulation with the Pb wall. 

K^2 decay followed by /i ^ e decay: energetic forward electrons contributing to background 
are suppressed according to the Michel distribution as muons from K^2 decays are fully polarised. 
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Figure 2: (a) Reconstructed squared missing mass distribution M^jj,j,(e) for the candidates: data (dots) 
presented as sum of signal and backgrounds (filled areas), (b) Numbers of Kf.2 candidates and background 
events in track momentum bin. 



Source 




Source 


NB/N,or 


Source 


Ns/Ntot 




(6.28 ±0.17)% 


Kely (SD) 


(1.02 ±0.15)% 


Ke3 


0.03% 


K^t2 (M ^ e) 


(0.23 ±0.01)% 


Beam halo 


(1.45 ±0.04)% 


Km 


0.03% 



Total background: (8.03 ± 0.23)% 



Table 1: Summary of the background sources in the Kgi sample. 

Ke2y (SD) decay, a background by Rk definition, has a rate similar to that of Ke2'- experimen- 
tally BR = (1.52±0.23) X 10^^ [^. Theoretical BR calculations have similar precision, depending 
on the form factor kinematic dependence model. Energetic electrons {E* > 230 MeV in K frame) 
with 7 escaping detector acceptance contribute to the background. MC background estimation has 
a 15% uncertainty due to limited knowledge [^] of the process. A recent precise measurement by 
KLOE iQ], published after announcement of the NA62 preliminary result, is not taken into account. 

Beam halo background in the Ke2 sample induced by halo muons undergoing jx e decays in 
flight is measured directly using the K data samples. Background rate and kinematical distribution 
are qualitatively reproduced by a muon halo simulation. The uncertainty is due to the limited size of 
the K sample. Beam halo is the only significant background source in the K^2 sample, measured 
to be 0.25% (with a negligible uncertainty) with the same technique as for Ke2 decays. 

The number of K(2 candidates is N{Ke2) = 51,089 (about four times the statistics reported by 
KLOE [|]) and (^^2) = 15.56 X 10^. The M^jj,j,(e) distributions of data events and backgrounds 
are presented in Fig. Backgrounds integrated over track momentum are summarized in Table [I]; 
their distributions over track momentum are presented in Fig. 

4. Systematic uncertainties and results 

Electron identification efficiency is measured directly as a function of track momentum and 
LKr impact point using pure samples of electrons obtained by kinematic selection of — > 7r"e+v 
decays collected concurrently with the Ke2 sample, and Kf^ —)■ n^e^v decays from a special 15h 
run. The and measurements are in good agreement. The measured fg averaged over the Ke2 
sample is 99.2%, with a precision better than 0.05%. Muon identification inefficiency is negligible. 
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Figure 3: (a) Rk in track momentum bins, (b) Summary of Rk measurements. 



Source 


5Rk X 10^ 


Source 


5Rk X 10^ 


Source 


5Rk X 10^ 


Statistical 


0.012 


Beam halo 


0.001 


Geom. acceptance 


0.002 




0.004 


Electron ID 


0.001 


Trigger dead time 


0.007 


Kely (SD) 


0.004 


IB simulation 


0.007 







Table 2: Summary of uncertainties of R^'- statistical and systematic contributions. 

The geometric acceptance correction A{K^2) / ^{Kei) depends on the radiative K^jj (IB) 
decays. A conservative systematic uncertainty is attributed to approximations used in the Kf.2y IB 
simulation, which follows [|l|]. The resummation of leading logarithms is not taken into account, 
however no systematic error is ascribed due to that. An additional systematic uncertainty reflects 
the precision of beam line and apparatus description in the MC simulation. 

Trigger efficiency correction e{Ke2) / £{K^2) ~ 99.9% accounts for the fact that Ke2 and K^2 
decay modes are collected with different trigger conditions: the £" > 10 GeV LKr energy deposition 
signal enters the Ke2 trigger only. A conservative systematic uncertainty of 0.3% is ascribed due to 
effects of trigger dead time which affect the two modes differencly. LKr global readout efficiency 
/lkt is measured directly to be (99.80 ±0.01)% using an independent LKr readout. 

The independent measurements of Rk in track momentum bins with are presented in Fig. ^a. 
The uncertainties are summarised in Table |[ The preliminary NA62 result is Rk = (2.500 =b 
0.012stat. iO.Ollsyst.) X 10"^ = (2.500 ±0.016) x 10"^ consistent with the SM expectation. The 
whole 2007-08 data sample will allow pushing the uncertainty of Rk down to 0.4%. A summary 
of Rk measurements is presented in Fig. the new world average is (2.498 ±0.014) x 10^^. 
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